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ABSTRACT 

Motivated by the suggestion of Kang, Ryu & Jones (1996) that particles can be acceler- 
ated to high energies via diffusive shock acceleration process at the accretion shocks formed 
by the infalling flow toward the clusters of galaxies, we have calculated the expected particle 
flux from a cosmological ensemble of clusters. We use the observed temperature distribution 
of local clusters and assume a simple power-law evolutionary model for the comoving den- 
sity of the clusters. The shock parameters such as the shock radius and velocity are deduced 
from the ICM temperature using the self-similar solutions for secondary infall onto the clus- 
ters. The magnetic field strength is assume to be in equipartition with the postshock thermal 
energy behind the accretion shock. We also assume that the injected energy spectrum is a 
power-law with the exponential cutoff at the maximum energy which is calculated from the 
condition that the energy gain rate for diffusive shock acceleration is balanced by the loss 
rate due to the interactions with the cosmic background radiation. In contrast to the earlier 
paper we have adopted here the description of the cosmic ray diffusion by Jokipii (1987) 
which leads to considerably higher particle energies. Finally the injected particle spectrum at 
the clusters is integrated over the cosmological distance to earth by considering the energy 
loss due to the interactions with the cosmic background radiation. Our calculations show that 
the expected spectrum of high-energy protons from the cosmological ensemble of the cluster 
accretion shocks could match well the observed cosmic ray spectrum near 10 19 eV with rea- 
sonable parameters and models if about 10 -4 of the infalling kinetic energy can be injected 
into the intergalactic space as the high energy particles. 
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1 INTRODUCTION 

It is widely believed that diffusive shock acceleration is the mech- 
anism from which cosmic rays get their energy. There are vari- 
ous models to account for the origin of cosmic rays below about 
3xl0 18 eV = 3EeV. At low energies, up to about 10 14 eV, su- 
pernova explosions into the interstellar medium give a reasonable 
and successful explanation for the data for protons. At higher en- 
ergies, several models have been proposed, most notably a galac- 
tic wind termination shock (Jokipii & Morfill 1987), and multiple 
shocks in an ensemble of OB superbubbles and young supernova 
remnants (Axford 1992). A comprehensive, albeit tentative, the- 
ory has been proposed by Biermann (1993 and later papers) that 
explains the cosmic ray spectrum, with its chemical abundances 
and the knee feature as resulting from a combination of supernova 
explosions into the interstellar medium, and supernova explosions 



into strong stellar winds of progenitor Wolf-Rayet stars (Biermann 
1993; Biermann & Cassinelli 1993; Biermann 1995). Above the 
"ankle" at about 3 EeV (the ultra-high energy cosmic ray regime, 
called UHECR hereafter), a simultaneous change in spectrum and 
composition of the cosmic ray spectrum (Bird et al. 1994) suggests 
a change of origin. As argued already by Cocconi (1956) particles 
of higher energy need to come from outside our galaxy due to the 
very large Larmor radius, and therefore need to be accelerated in 
extragalactic sources. Here we concentrate on these high energy 
particles that almost certainly come from extragalactic sources. 

It is well known that the extragalactic cosmic ray spectrum 
must show the Greisen-Zatsepin-Kuzmin (GZK) cutoff at about 
60 EeV due to interactions with the cosmic background radiation 
(CBR), regardless whether protons, heavy nuclei or photons are 
considered as the energetic particles (Greisen 1966; Zatsepin & 
Kuzmin 1966). Pair production and the cosmological evolution of 
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radiation backgrounds set limits for the distance of cosmic rays 
even at lower energies, thus any extragalactic model of cosmic 
ray origin must propose at least some sources in our cosmologi- 
cal neighborhood which can account for the highest energies. In 
an expansion of an earlier proposition by Biermann & Strittmatter 
(1987) which predicted maximum proton energies in active galax- 
ies of 10 21 eV, Rachen & Biermann (1993) developed a model to 
accelerate cosmic rays up to a few 100 EeV at strong shocks at the 
end of extended jets in powerful radio galaxies (FR-II radio galax- 
ies, Fanaroff & Riley 1974). It has been shown that even for a mod- 
erate proton content in the jets this model can generally account 
for the cosmic ray flux above the ankle, and is consistent with air 
shower data suggesting a takeover from heavy to light nuclei in this 
energy range (Rachen, Stanev & Biermann 1993). It can be shown 
that this model can be extended to a maximum particle energy in 
the source of about 4xl0 21 eV (Biermann 1996). One prediction 
of this particular model has been tested successfully, and that is the 
expected correlation of arrival directions of very high energy cos- 
mic rays with the large scale distribution of radio galaxies, in the 
supergalactic plane (Stanev et al. 1995; Hayashida et al. 1996). 

However, the uncertain proton content of the jets, the energy 
limitation for cosmic rays due to the relatively small acceleration 
region in hot spots, the finite life time of hot spots, and the large 
distance of the closest well known FR-II radio galaxy set strong 
constraints on the predictability of the model, and leaves room for 
other contributions to the highest energy cosmic rays. In particular, 
after the detection of a 320 EeV air shower by Fly's Eye (Bird et 
al. 1994), and another 200 EeV event by the AGASA ground ar- 
ray (Hayashida et al. 1994), the various acceleration models have 
been critically reviewed (e.g. Elbert & Sommers 1995; Biermann 
1996) and new models have been suggested. Most of them are 
based on astrophysical objects whose physical properties are under 
controversial discussion, as the decay of topological defects (Bhat- 
tacharjee 1991; Sigl, Schramm & Bhattacharjee 1994; Protheroe & 
Johnson 1996; Protheroe & Stanev 1996), or rapid acceleration in 
Gamma Ray Burst sources (Milgrom & Usov 1995; Vietri 1995; 
Waxman 1995), but the poor statistics of events above 100 EeV al- 
lows a variety of other explanations, including the earlier proposal 
of radio galaxy origin (Rachen 1995; Stanev et al. 1995; Biermann 
1995, 1996). 

Recent work shows that there could be sites for shock ac- 
celeration in extragalactic space alternative to the strongly con- 
strained acceleration in radio galaxy hot spots. According to hydro- 
dynamic simulations of large scale structure formation (e.g. Kang 
et al. 1994a; Cen & Ostriker 1994), accretion shocks are formed 
in the baryonic component around non-linear structures collapsed 
from the primordial density inhomogeneities as a result of gravita- 
tional instability. Those structures can be identified as pancake-like 
supergalactic planes, still denser filaments, and clusters of galax- 
ies which form at intersections of pancakes, in any variants of the 
many cosmological models. They are surrounded by the hot gas 
heated by the accretion shocks and the particles can be accelerated 
to very high energies at these shocks via first order Fermi process. 
Kang, Jones & Ryu (1995) and Kang, Ryu & Jones (1996, KRJ96 
hereafter) suggested that the accretion shocks around the clusters 
of galaxies could be as fast as 1000-3000 kms _1 and so could be 
good acceleration sites for the UHECRs up to several 10 EeV, pro- 
vided there is a turbulent magnetic field so that the diffusion is in 
the Bohm limit, and if the magnetic field around the clusters is or- 
der of microgauss. We note here the maximum energy could be 
shifted to super-GZK values, if the field geometry is close to being 
quasi-perpendicular (Jokipii 1987). The significance of the accre- 



tion shocks is that they are the largest and longest lived shocks in 
the universe, so that they naturally get past the primary factors nor- 
mally limiting Fermi acceleration to such high energies. Especially 
the accretion shocks formed around the clusters of galaxies which 
have the deepest gravitational potential well are the strongest and 
thus could accelerate the particles to the maximum possible energy. 
Independently Norman, Melrose & Achterberg (1995) showed with 
the help of a similar argument that the shocks associated with the 
large-scale structure formation could accelerate the protons up to 
= 50 EeV if there is a primordial field of 1-10 nanogauss, 
or if microgauss field can be self-generated in shocks. 

In the present study we have estimated the contribution of the 
CR protons from an ensemble of the cluster accretion shocks dis- 
tributed in the universe by adopting some simple models for the 
accretion flows onto clusters, the strength, geometry, degree of ir- 
regularities of the magnetic field near the cluster accretion shocks, 
and the cosmological evolution of the cluster distribution. The de- 
tails of the models will be given in §2. In §3 the estimated CR pro- 
ton spectrum has been compared with actual observations in order 
to see if the proposed origin could explain the energy spectrum of 
UHECRs with reasonable physical parameters and models. We dis- 
cuss the implication of these results and models also in §3. 

We adopt for the following Qo = 1 and write the Hubble 
constant as Ho = 100 h km/s/Mpc; in numerical calculations we 
use h = 0.75. 



2 MODELS 

2.1 Accretion Shocks around Large Scale Structures 

It is generally accepted that both galaxy distribution in the ob- 
served Universe and matter distribution in numerically simulated 
universe in most of generic cosmological scenarios show sheet- 
like and filamentary structures on large scale (e.g. de Lapparent, 
Geller & Huchra 1991; Melott 1987; White et al. 1987). Numeri- 
cal simulations based on even a hierarchical clustering model such 
as variants of CDM models also indicate that the dominant non- 
linear structure is a network of filaments (White 1996) and the fil- 
aments became longer, straighter and clumpier as the hierarchical 
collapse proceeds (Summers 1996). Clusters of galaxies, especially 
the richest ones, form mostly at the vertices where several filaments 
intersect. They grow nonlinearly by gravitationally attracting mat- 
ter along the filaments and also strengthen and align the filaments at 
the same time. The dominance of filaments in the visual impression 
of large scale matter distribution could be understood in part by the 
facts that the filaments as mostly intersecting curves between pan- 
cakes should naturally have higher density contrast than pancakes 
themselves (see West, Villumsen & Dekel 1991), and that a simu- 
lated universe is most often presented as a two-dimensional projec- 
tion or a slice cut where sheetlike structures cannot be easily seen. 
It also depends to some degree on choosing a right level of density 
contrast. Thus, for example, knots (clusters) will become dominant 
if one chooses to look at only the highest density peaks. 

Although formation of shocks and subsequent heating of the 
gas by them, when matter accretes toward these large-scale co- 
herent structures, are implicit in any cosmological hydrodynamic 
simulations, the visual impression of accretion shocks has became 
clear only after high resolution Eulerian codes were introduced 
to the numerical cosmology (e.g. Ryu et al. 1993; Bryan et al. 
1995). Interested readers are referred to Kang et al. (1994b) for 
the details of a comparison of various cosmological hydrodynamic 
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Figure 1. A representative slice cut of a simulated universe in a standard cold dark matter model. The left panel shows the X-ray luminosity distribution while 
the right panels shows the gas temperature distribution. See text for the details. 



codes. They have demonstrated the accretion shocks can be cap- 
tured clearly in the simulations done by these high-resolution Eu- 
lerian codes, while their existence is not so obvious in the simu- 
lations done by the particle codes based on the SPH method (see 
their Figs. 5). Fig. 1 shows a slice cut through a simulated universe 
based on the standard CDM model previously reported by Kang 
et al. (1994a). Readers are referred to their paper for the details of 
the cosmological model parameters, which should not be crucial to 
the current discussion. The left panel shows the X-ray luminosity 
modeled by L x = p^T 1 ^ 2 and could be a representative distri- 
bution of X-ray clusters. Being weighted toward the higher den- 
sity region, this shows the knot-like distribution with some align- 
ments of several knots into filaments/sheets. The right panel shows 
the gas temperature distribution. The hotter region of T > 10 6 K 
is represented by solid contour curves, while the colder region of 
T < 10 6 K by dotted contour curves. Shocks can be seen clearly 
as strong gradients in temperature distribution, and they encompass 
the region of moderately overdense regions of sheetlike/filamentary 
structures. The most prominent cluster in the lower-left corner has 
been shown in Fig. 1 of KRJ96 in which the velocity pattern clearly 
shows accretion flows from the background IGM and along the fila- 
ments toward the cluster. These figures demonstrate that the accre- 
tion shocks associated with the large-scale coherent structures do 
indeed exist, even though any direct observations of such shocks 
have not been made yet. It is apparent from the diffusive shock ac- 
celeration theory that some particles are accelerated by these large- 
scale shocks. The accretion shocks around clusters of galaxies, be- 
ing the fastest shocks generated by the deepest potential wells, are 
our focus in the present paper. 



2.2 Self-Similar Evolution of Clusters 

The temperature of the intracluster medium (ICM) within a clus- 
ter can be obtained with a reasonable accuracy from the X-ray ob- 
servations (e.g. David et al. 1993). Fortunately, important physical 
parameters such as the velocity of an accretion shock, the veloc- 
ity dispersion of galaxies, and the depth of the potential well of 
the cluster can be related rather well with the ICM temperature. 
Thus one can deduce such parameters from the observed temper- 
ature of X-ray clusters, for example, by assuming that the gas is 
in a hydrostatic equilibrium with the gravitational potential due 
to total mass including both the dark and baryonic matter. Here 
we take another approach in which one-dimensional (ID) spheri- 
cal collapse of clusters is modeled as the collapse of an initially 
overdense point-mass perturbation followed by the secondary in- 
fall of background medium. Such models have been studied both 
semi-analytically (Fillmore & Goldreich 1984; Bertschinger 1985) 
and numerically (Ryu & Kang 1996). In such an accretion flow the 
infalling baryonic matter is stopped by a shock, while the collision- 
less matter forms many caustics. In this paper we assume that the 
evolution of clusters and the properties of the accretion shocks can 
be represented by such flows. Although this ID approach does not 
account for the virialization of the central region of clusters, we will 
show below essential characteristics of real clusters can be related 
with parameters of the accretion shock through this model. 

One can expect the evolution of such a flow and its accre- 
tion rate should be dependent upon the expansion of background 
universe for the given initial power spectrum of density pertur- 
bations (i.e. P(k) is constant for a point-mass perturbation). The 
flow solution in an f2 = 1 universe approaches a self-similar 
form (Bertschinger 1985) due to the scale-free nature of the prob- 
lem and so can be treated analytically. In a low density universe 
(i.e. fi < 1), however, either open or flat with the cosmological 
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constant, the solution is not self-similar and can be studied only 
numerically (Ryu & Kang 1996). Here we will consider the ac- 
cretion flows only in the Einstein-de Sitter universe (Qo = 1), 
and adopt the one-dimensional (ID), self-similar accretion given 
by Bertschinger (1985) as an evolutionary model for clusters and 
the accretion shocks. 

Since the solution is self-similar, the shock parameters such 
as the radius and velocity of the accretion shock, and the post- 
shock gas temperature can be uniquely determined by a single 
parameter at a given epoch. For example, the mass contained in- 
side the outermost caustic, M c = M(r < r c ) at present epoch 
can be such a parameter. The radius and velocity of the accretion 
shock, and the postshock temperature at r = 0.3r s = 0.64/i -1 
Mpc are given by r s = 2.12/i _1 Mpc(M c h/10 15 M ) 1/3 , 
V s = 1.75xl0 3 kms _1 (M c h/10 15 M G ) 1/3 , and T c (0.3r s ) = 
6.06keV(M c ft/lO 15 M ) 2/3 , respectively (Ryu & Kang 1996). 
The reason that T c (0.3r s ) is an interesting quantity will be given 
shortly. If one follows the evolution of a given perturbation, on the 
other hand, the length scale of the accretion flow which is pro- 
portional to the turn-around radius will grow with time as r ta oc 
t 8 '" i oc (1 + 2)~ 4 / 3 . The position of the shock, r s , in units of r ta 
and the velocity of the shock, V s , in units of r ta /£ are fixed, so 
the shock radius and velocity for the same perturbation evolves as 
r s oc (1 + z)- A ' 3 and V s oc (1 + z) 1 ' 6 . 

Here we attempt to relate the physical parameters of the accre- 
tion shocks with observed temperature of X-ray clusters which of- 
ten represents the emission weighted temperature in the core within 
~ 0.5/i _1 Mpc. We note that the temperature of the postshock gas 
in the ID self-similar solution increases toward the center, in fact, 
to infinity, while observed ICM temperature distribution of most 
clusters is nearly isothermal, and in some cases shows a clear de- 
pression towards the center due to cooling. This discrepancy comes 
about, because the extrapolation of the self-similar evolution to 
t — > U (r — > 0) is not valid during the initial collapse when the 
accreted mass is not much greater than the mass inside the initial 
perturbation. Also the ID self-similar solution does not account for 
the virialization of the core region. According to the 3D cosmo- 
logical hydrodynamic simulations (e.g. Crone, Evrard & Richstone 
1994; Kang et al. 1994a; Navarro, Frenk & White 1995), the ICM 
is shock heated to the virial temperature and then settles into hydro- 
static equilibrium with an approximate isothermal structure, which 
is consistent with the observed temperature distribution of the clus- 
ters. Cosmological SPH simulations (Navarro et al. 1995; Evrard, 
Metzler & Navarro 1995) showed that simulated clusters of differ- 
ent masses in fact have similar structures when scaled to a fixed 
density contrast (e.g. 5(r) = p(r)/p cr it ^ 200—500). Navarro et 
al. also showed that the temperature profile of simulated clusters 
can be approximated within a factor of two by that of the ID self- 
similar solutions for the outer region, r J> 0.3r2oo (where V2oo is 
the radius at 5 — 200), while the inner region, r < 0.3r2oo, is 
isothermal. From this consideration, we have made a simple ap- 
proximation that the observed X-ray temperature of a cluster is 
similar to T c (0.3r s ) of the self-similar solutions, and so it provides 
the direct information about the shock velocity and radius from the 
self-similar solutions. 

From the equations given above which relates r s , T c (0.3r s ), 
and V s with the cluster mass M c and their redshift dependences, 
one can find them as a function of kT oba = T c (0.3r s ) and z in an 
fio = 1 universe as follows. 

*■■ = 2 ' m " MpC (6^) 1/2 ( 1 + ^ 3/2 (1) 



V e = 1.75X10- km.- (^) . (2) 

Evrard et al. , on the other hand, showed that rsoo where <5(r5oo) = 
500 is a conservative estimate for the boundary separating the inner 
virialized region from the outer Mailing flow. They gave the scal- 
ing relation r 5 oo = 0.965 h' 1 Mpc x (fcT/6.06 koV) 1/2 for clus- 
ters at z = for an Qq = 1 universe. The fact that the ratio of the 
shock radius of the self-similar solution to the characteristic radius 
of virialized region in 3D simulation is r s /r5oo = 2.2 seems rea- 
sonable, since r s corresponds to the radial position at S ~ 80 and 
so it should be larger than rsoo- According to Ryu & Kang (1996), 
the clusters of a given temperature have smaller accretion veloc- 
ity by less 20% for Qo = 0.3 in either open or flat with Ao 7^ 
universes. Since we use observed temperature distribution for the 
cluster abundance to be discussed below, the main results of our 
model should remain valid for a low flo cosmology. 

According to Bertschinger (1985), the gas density upstream 
to the shock is pi = 4.02fi b Pait(z) = 7.56x 10~ 29 h 2 g cm" 3 x 
fib(l + z) ' i . Thus all the shock parameters necessary for our model 
can be obtained from the observed redshift and X-ray temperature. 
In the following we adopt VL b = 0.06 in all numerical calculations. 



2.3 Magnetic Field and Diffusion Coefficient 

The strength and morphology of the intergalactic magnetic fields 
remain largely unknown and observational tasks to detect them are 
very challenging even with today's technology (see, for a review, 
Kronberg 1994). Theoretical study on the generation of the primor- 
dial fields and their subsequent amplification during the structure 
formation is also still in its infancy. In fact it is not clear at all that 
there was a primordial magnetic field. Recent study by Kulsrud et 
al. (1996), however, showed that a weak seed field can be generated 
at shocks and then amplified via the protogalactic turbulence during 
the structure formation up to equipartition with the turbulence. An 
alternative theory (Biermann 1996) uses magnetic stellar winds as 
the sources of magnetic fields for galaxies and their environment. 
On observational fronts, there exist some concrete observations that 
can be used in inferring the general distribution of the magnetic 
fields in intergalactic space. The current observational upper limit 
for a large-scale, pervading field is about 10~ 9 r gauss (Kron- 
berg 1994), where ro is the assumed reversal scale of the magnetic 
field topology in units of 1 Mpc. Thus, if the reversal scale through- 
out the universe were larger than 1 Mpc, such as « 30frT 1 Mpc, the 
bubble-scale of the galaxy distribution, then this upper limit would 
be reduced by a factor of 6.3 (for h = 0.75). The magnetic fields 
inside typical galaxies are observed to be order of 3-10 /igauss, and 
magnetic fields near the /igauss level are also common in core re- 
gions of rich clusters according to many recent observations (Kim, 
Tribble & Kronberg 1991; Taylor & Perley 1993; Taylor, Barton & 
Ge 1994); some cooling flow clusters clearly show magnetic fields 
even higher than those typical in the interstellar medium of galax- 
ies (see Kronberg 1994). Studies on the field generation via dy- 
namos in the cooling flows (Ruzmaikin et al. 1989) and the field 
inputs from radio galaxies inside clusters (EnBlin et al. 1996) have 
attracted some attention recently. Here we are concerned most with 
the fields near the accretion shocks on larger scale than cluster 
core, that is, within ~5/i _1 Mpc around clusters of galaxies. Ob- 
servation of such fields has been attempted by Kim et al. (1989) in 
which fields of 0.1 /igauss level were deduced in the plane of the 
supercluster connecting the Coma cluster and A1367 by assuming 
equipartition between the magnetic fields and relativistic particles. 
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This "bridge" of emission region of ~ 1.125/i Mpc seems to be 
infalling toward the Coma and is within the supergalactic plane (so 
presumably inside the accretion shocks). Vallee (1990,1993), on 
the other hand, suggested the existence of a magnetic field compo- 
nent of ~ 1.5 /igauss in a 10 Mpc region around the Virgo cluster, 
in this case using an assumed length scale of 10 Mpc. This might 
represent the fields outside the accretion shocks encompassing the 
supergalactic plane. 

Here we adopt a simple, but common assumption that the ICM 
magnetic field is in rough equipartition with the thermal energy of 
the gas. Then the field strength in postshock region can be estimated 
from the thermal energy inside the shock (i.e. in the postshock re- 
gion) according to 



1.5 



Pi 
pm H 



kT 2 



8tt 



(3) 



where p 2 (z) = 4pi(z), p, = 0.61, and T 2 = 2.69xl0 7 K x 
(fcT o bs/6.06keV). Then the postshock field is given by 



B 2 — (1.71 /igauss) fsh 



k-T 



0.06 6.06keV 



(4) 



where fs Si 1 is a factor which controls the field strength in terms 
of the equipartition value, that is, f B — 1 means the equiparti- 
tion between the field and thermal energies. For fiducial values of 
n b = 0.06 and h = 0.75, this gives for f B = 1 the fields of 0.52- 
1.64 /igauss for fcT b s = 1— lOkeV at z = 0. It is much harder 
to estimate the fields in the unshocked infalling flow upstream to 
the shock, so here we will simply assume B\ = Bz(p\/p<2) for a 
turbulent field. We note here that the turbulent amplification of the 
field could generate some fields not only downstream but also up- 
stream to the shocks due to turbulent nature of the flows (Kulsrud 
et al. 1996). 

We consider two kinds of models for the particle diffusion. 
The theoretical minimum for the diffusion coefficient in a strongly 
turbulent field with parallel geometry (i.e. the mean field is parallel 
to the flow direction) is given by the Bohm formula (see, e.g. , 
Drury 1983). 



KB 



3 ' 



(5) 



where r g is the gyroradius and v is the velocity of the particle. 
On the other hand, the minimum diffusion coefficient in the per- 
pendicular shocks (i.e. the mean field is perpendicular to the flow 
direction) (Jokipii 1987) is given by 



kj = r g V B = 3 (Vb/v) KB' 



(6) 



For high energy particles, v ~ c, so kj is smaller than kb by a 
factor of Vs/c. 

The rate of diffusive acceleration is determined by the diffu- 
sion coefficients k = K\\ cos 2 + k± sin 2 0, where K\\ and k± 
are the diffusion coefficient parallel and perpendicular to the mag- 
netic field, respectively, and is the angle between the field and the 
shock normal (Jokipii 1987). Here we will follow Jokipii (1987) 
in assuming that k±/k\\ = (1 + ?/ 2 ) -1 and K\\ = i]r s v/3, where 
r\ is the ratio of the mean free path parallel to the magnetic field 
to the gyroradius. This is the result from standard kinetic theory 
which does not take account of field-line meandering. The Bohm 
diffusion coefficient corresponds to the case where the field is tur- 
bulent on a scale of r g at all momenta, so it is effectively equiva- 
lent to the case of r\ 1 and n± ~ k\\ with a random distribution 
of the obliquity. On the other hand, the minimum diffusion coef- 
ficient in a perpendicular shock, which is referred as the Jokipii 



diffusion throughout this paper, corresponds to the case of 8 = 90° 
and r/ max ~ c/(SV s ) = 100/14,3, where V s ,3 is the shock veloc- 
ity in units of 10 3 km s _1 . The latter is derived from the condition 
that Kx jm i n ~ fgVs ~ (r/r g c)/[3(l + rf)\, which implies that the 
particles should be scattered before they drift through the shock. 
This condition is necessary to maintain the isotropy of the particle 
distribution. 



2.4 Energy Losses and Maximum Energy 

The mean acceleration time scale for a particle to reach a momen- 
tum p is determined by the velocity jump at the shock, and the 
diffusion coefficient (e.g. , Drury 1983; Jokipii 1987), that is 



X 



X ~ 1 V. Vs 
where \ is the compression ratio of the shock and 



(7) 



2 , sin 
cos + 



1 + 7] 2 



' + x 2 sin 2 0/(1 + if) 



[cos 2 + x 2 sin 2 6<] 



3/2 



In the strong shock limit (i.e. \ = 4), it can be written as 



(4.23x10 years) 



B»V 2 - 



.7(4,0), 



(8) 



where E\s is the particle energy in units of EeV, B M is the field 
strength in units of microgauss. 

The protons lose energy due to pair production and pion pro- 
duction on the CBR not only on their way to earth, but also during 
their acceleration at the cluster shocks. At energies around 10 EeV 
the energy loss due to pair production is dominant and the loss time 
scale can be as short as (5x 10 9 years)/(l + z) 3 . Above the GZK 
cutoff at about 60 EeV photopion production becomes dominant 
and the loss time scale can be as short as (5x 10 7 years) /(l + z) 3 . 
The maximum energy -E max up to which the protons can be acceler- 
ated by the cluster accretion shocks is found by setting r acc = Tint, 
with ri n t being the time scale for interaction losses with the CBR, 
given by 



Tint = 



C 7 



kQ(z) 

de e (<7k) (e) In 



(9) 



exp ( 



2-ykQ(z) 



where 7 is the Lorentz factor of the cosmic ray protons, Q(z) = 
0o (1 + z) the temperature of the cosmic microwave background at 
epoch 2, and (o~k) the inelasticity weighted cross section of the pj 
reaction, averaged over all final states (Rachen & Biermann, 1993). 
Equation (Q) considers both pair production and pion production 
on the CBR; the impact of a putative cosmic infrared background 
is negligible for the determination of i? ma x. 

In Fig. 2 we have shown the loss time scales due to the in- 
teractions with CBR at z = 0, and the acceleration time scales 
in both Bohm and Jokipii diffusion limits. Here the acceleration 
time scales for V B = 1000kms~ 1 and Bi = 1/xgauss are plotted. 
For other values of V s and B, the acceleration time scales can be 
found by the following scaling relations: r acc ,B oc V S ~ 2 B~ 1 for the 
Bohm limit, and r acCj j oc V S ~ 1 B~ 1 for the Jokipii limit. The max- 
imum energy i? max for a shock of these parameters can be found 
by the intersection of two curves of Ti nt and r acc . The characteris- 
tic shape of loss time scale Tint as a function of energy causes only 
a weak dependence of i5 max on r acc , if r acc <^ 3x 10 9 years, but 
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Figure 2. Schematic representation of acceleration and energy loss time 
scales as functions of particle energy. The dotted and dashed lines repre- 
sent the loss time scale due to pion-production and pair-production, respec- 
tively, on both CMB and infrared background radiation. The heavy solid 
line is the added loss time scale due to both interactions, and it turns out 
that the infrared contribution is nowhere relevant. The dot-dashed and light 
solid line represent the acceleration time scales for Bohm and Jokipii dif- 
fusion models, respectively. For this case, adopted shock parameters are 
Vs = 1000 km s _1 and B = 1/^gauss. The intersection point of both curves 
determines the maximum energy of acceleration. 



a strong dependence if the acceleration is slower. This time scale 
corresponds to n a t for E ~ 10 19 ' 6 eV, Thus one can see from Fig. 
2 that, for a canonical accretion shock, the maximum energy for the 
Bohm limit is most likely to be set by the pair-production loss and 
ranges from 1 to lOEeV, while that for the Jokipii limit is set by 
the pion-production loss at around 50 EeV. 

According to equation the acceleration time scale is in 
general dependent upon the obliquity and the strength of turbulent 
field. Fig. 3 shows the acceleration time scales for the particles to 
reach up to E = 10 19 8 eV for r) = 1, 10 and 100 as a function 
of the upstream oblique angle 6. Once again the acceleration time 
scales for V s = 1000 km s~ and for B\ = l/_tgauss are plotted. 
This shows that the particles can be accelerated above 10 19 ' 6 eV 
only when the mean field is nearly perpendicular to the shock nor- 
mal. The obliquity dependence of r acc is rather weak in the strong 
scattering limit (77 ~ 1), while it becomes very strong when the 
cross-field diffusion is small (77 ~ 100). One can also note that r aC c 
linearly increases with 77 in quasi-parallel shocks (i.e. 6 m 0°), so 
the acceleration is most efficient in the limit of strong turbulences 
(i.e. 77 — > 1). On the other hand, r acc is nearly inversely propor- 
tional to 77 in quasi-perpendicular shocks (i.e. 6 w 90°), so the 
acceleration becomes most efficient in the limit of weak scattering 

(i.e. 77 — > ?7max). 



Figure 3. Time scales for the protons to be accelerated to E = 10 19 ' 6 eV 
by a shock of V B = 1000kms _1 and Bi = 1/j gauss as a function for 
obliquity. The curves are labeled with the values of r\ = A/r g = 1, 10, 
and 100. 



2.5 Injection Spectrum 

For each cluster with given values of fcT b s and 2, the self-similar 
solution gives the shock parameters, pi, V s , r s , and B\, and the 
time scale condition in equation (^) gives the estimate for the cut- 
off energy in the particle spectrum. Momentum distribution of the 
protons at the shock is assumed to be a power law whose index 
is given by the parameter a, and with an exponential cutoff at the 
maximum energy (i? m ax = cp c ). The minimum momentum for 
the power law is pi n j ~ m p V s . The power-law index is a — 4 in 
the limit of strong shocks with dynamically insignificant CR en- 
ergy density (i.e. test-particle limit), but can vary slightly around 
this value for real shock waves (we note that a — 4 corresponds 
to a power-law in the energy spectrum of E~ 2 ). Here we consider 
values of a equal to and slightly larger than 4. We further assume 
that the process with which the particles escape from the shock is 
momentum-independent, that is, the particle spectra injected into 
intergalactic (IG) space has the same shape of the proton spectra at 
the shock. Then the injected distribution from a cluster with a given 
ICM temperature is given by 



f(T,z,p) = A(T,z)p- a e X p 



V 



MT,z) 



(10) 



where A(T, 2) is a normalization constant. It is chosen by assuming 
that a small fraction (e) of the kinetic energy density of the infalling 
matter in the shock rest frame (piV 2 ) is converted to CR energy 
and then injected into IG space. So the free parameter, e, along with 
the assumed value of a controls the amplitude A(T, z) according 
to 



4ty 



f(T,z,p)Ep 2 dp 



(11) 



Pinj 
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Thus, for given values of e and a, the normalization constants were 
numerically calculated for given values of T and z. For a = 4, this 
gives an approximation for A(T, z) such as 



A(T,z) 



epiV s 2 



4:TTm p c 2 In (pc/mpc) 



(12) 



This along with p c = (1/c) E max (T, z) completely specifies the 
injection spectrum from a cluster. 

In KRJ96, it was assumed that the particles escape from the 
shock once they diffuse upstream to a distance comparable to the 
radius of the shock due to the lateral diffusion. Since this escape 
process depends on the diffusion length, the particle spectrum in- 
jected into the IG space has a strong dependence on the momentum, 
that is, only highest energy particles with longest diffusion length 
can diffuse out from the shock. We note that, if we take the same 
momentum-dependent escape model, the resulting spectrum cannot 
fit the observed CR distribution over a wide range of the particle en- 
ergy as well as the spectrum of our momentum-independent escape 
model can. We will also discuss the escape process in §3.1. 

2.6 Cluster Distribution Function 

The evolution of the cluster population is a critical issue which 
has been under intense discussion recently, since it can provide a 
potentially powerful tool for discriminating different cosmological 
models, especially the matter density of the universe (i.e. £7q). The- 
oretical prediction based on a hierarchical clustering model is that 
clusters are fainter but more abundant in the past (Kaiser 1986). The 
evolution of cluster population can be calculated by using various 
methods based on so-called Press-Schechter formalism, if the back- 
ground cosmology (i.e. fio and Ao) and the initial density power 
spectrum are given (e.g. Eke, Cole & Frenk 1996; Kitayama & 
Suto 1996; Bond & Myers 1996). The density fluctuations of clus- 
ter mass (M ~ 10 1 M@) would form earlier in low fio universe 
than in high flo universe, while their formation rate at the present 
epoch would be higher in low bias models (i.e. larger as, Cen & 
Ostriker 1994; Kitayama & Suto 1996; Bond & Myers 1996). The- 
ses studies as well as most numerical studies (Kang et al. 1994a; 
Tsai & Buote 1996) indicate that the standard CDM model of a 
critical density universe produces too many clusters compared to 
observed local cluster abundance. Also recent observations of dis- 
tant X-ray clusters at z < 0.3. (Castander et al. 1995; Ebeling et al. 
1995) seem to indicate that clusters have evolved very little or not 
evolved at all in this redshift range, which is more consistent with 
low bias or low density models. 

In order to model the distribution of cosmological population 
of clusters, first we adopted a temperature distribution function of 
clusters given by Henry and Arnaud (1991) which was derived from 
the observed local clusters. The number density in unit comoving 
volume at the present is given by for kT =l-10keV, 



(k l \ 
7-T7 I 
keV/ 



(13) 



Secondly, we made a simple assumption that the comoving density 
of cluster evolves as a power-law of a scale factor, (1 + z), that is, 
n(kT, z) = n (kT)(l + z) m . The value of m at redshifts z < 0.3 
is most likely zero or slightly negative according to the observations 
mentioned above. According to the scaling law of Kaiser (1986) for 
a scale free power spectrum in a critical density universe, for CDM 
like density power spectrum of n e ff ~ — 1 at the cluster mass scale, 
the index m would be —0.7 for the above power-law temperature 
distribution given by equation (|l3|). Here we will consider a range 



of values for m (— 1 < m < +1). We take the minimum value of 
the redshift as that of the Virgo cluster, that is, z m i n = 0.0036. 
Most X-ray bright, rich clusters would form after z ~ 5 and in fact 
clusters of M = 10 Mq form around zt ~ 0.6 in an Qq = 1 
universe. From these considerations, we set the maximum redshift 
to be z max = 5. But in fact contributions from the epoch earlier 
than 2 ~ 1 are insignificant even for m — 1 modeled because of 
the spatial dilution (1/r 2 factor) and the interactions with CBR for 
distant sources. 

Then number of clusters between z and z + dz for a given 
temperature kT is given by 

fdVo 



dN{kT,z) =n (kT) (l + z) r 



dz 



dz 



(14) 



where dV c /dz = 4n(c/Ho) 3 [l - (1 + z)~ 1/2 } 2 /{l + z) s/2 for 
Qo — 1 and (dV c /dz) dz is the comoving volume of a shell de- 
fined between z and z + dz. Adding up the contribution from clus- 
ters from 2min to Zmax while including the interactions with CBR, 
and integrating over the temperature distribution, one can get the 
particle flux observed at earth according to 



J(E) 




(15) 



M(E)f(T,z,E) 



d c i(z) — (6 ft. Gpc) (1 + z — y/1 + z) is the luminosity distance 
of a cluster at redshift z, given in an Qo — 1 cosmology. M{E) is 
the modification factor that accounts for the interactions with CBR 
along the pathway to earth, calculated as defined in Rachen and 
Biermann (1993), using the the continuous energy loss approxima- 
tion introduced by Berezinsky and Grigor'eva (1988). It has been 
pointed out that the characteristic spikes occurring in M{E) close 
to the cutoff are exaggerated because of the continuous loss approx- 
imation (Yoshida & Teshima 1993; Protheroe & Johnson 1996). In 
the integration over redshift, however, those features are smoothed 
out, and it should only be noted that the cutoff is not quite as sharp 
as proposed by this method. 



3 RESULTS AND DISCUSSION 

Fig. 4 shows the maximum energy calculated according to equa- 
tion r acc = ri n t for both diffusion models for clusters with kT — 
2, 4, 6, and 10 keV at z = 0.0 as a function of the parameter 
/s . As found in KRJ96, the maximum energy cannot go above the 
GZK-cutoff for the Bohm diffusion coefficient. Thus if the gen- 
eral field direction is radial in the accretion flow, then the diffu- 
sive acceleration is too slow to accelerate the protons above the 
GZK-cutoff. We note from Fig. 3, however, that the transition from 
quasi-perpendicular Jokipii diffusion to Bohm diffusion does not 
lead to a linear decrease of -E max , so that the model can work also 
for moderately oblique shocks with r\ 1. 

Fig. 5 shows the expected proton spectrum from the cluster 
accretion shocks calculated for both diffusion models. The power- 
law index for the cosmological evolution is m — 0, the magnetic 
field strength parameter is fs = 1 for all models. Different values 
of the injection energy fraction e is assumed for each value of the 
spectral index a to obtain the better fit with the observation around 
lOEeV for Jokipii models and around 1 EeV for Bohm models, 
respectively. As expected, higher e is required for steeper spectra. 
For Jokipii diffusion, the models with 4.0 < a < 4.2 show good fits 
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Figure 4. Maximum energy of the protons accelerated by the accretion 
shocks associated with the clusters with ICM temperatures of kT = 
2, 4, 6, and 10 keV at z = as a function of the parameter fg (il^, = 0.06 
and h = 0.75). The top panel shows models with Jokipii diffusion, the bot- 
tom panel shows models with Bohm diffusion. 



to the data, while the models with Bohm diffusion produce too few 
particles above 1 EeV. 

The sensitivity to the cluster evolutionary model is shown in 
the top panel of Fig. 6. For all models here a — 4, e = 5x 10~ 5 , 
and fs = 1. The power-law index for the evolution with the 
redshift for the comoving density of clusters are m = —1,0, 
and +1 for the dotted, solid, and dashed lines, respectively. Since 
the nearby clusters contribute most, the resulting spectrum is not 
severely dependent on the evolutionary model at high redshifts. 
The relatively weak cosmological evolution of galaxy clusters com- 
pared to radio galaxies provides an even better fit to the light com- 
ponent data derived from the Fly's Eye air shower analysis (Rachen 
et al. 1993), but we point out that the errors are large here and that 
these data have not been confirmed by the AGASA collaboration. 

The bottom panel of Fig. 6 shows how the results depend on 
the magnetic field strength. For all models here a = 4, e = 
5xlCT 5 , and m = is used. The field strength relative to the 
equipartition strength is represented by fs = 0.1,0.5 and 1 for 
the dashed, dotted, and solid lines. One can see that cluster ac- 
cretion shock cannot produce enough particles above 10 EeV, if 
the field strength is much smaller than the equipartition value (i.e. 
/s ^ 0.1). The reduction of the field from the equipartition by a 
factor of up to two, however, still can give a reasonable fit to the 
observations. 

3.1 The escape process 

We need to consider why nature would choose an injection effi- 
ciency at a level of order 10~ 4 . It is believed that a supernova 
blast wave can transfer a few (1—10)% of the explosion energy 
into the cosmic ray component and the particles accelerated by the 
remnant shock are injected into the interstellar medium when the 
shock either decays into a sonic wave or breaks up due to various 
instabilities. It would be natural to assume that a similar fraction 
of the kinetic energy of the Mailing flow is likely to be transferred 
to the CRs. Unlike a supernova blast wave which expands from a 
point-like explosion into the interstellar space, however, an accre- 
tion shock forms due to the converging flow associated with grav- 
itational collapse. The particles will be mostly confined near the 
shock and so the escape of particles by swimming against infalling 
flows upstream will be very rare. This calls for some physical mod- 
els of escape of CRs from the cluster shocks which can explain that 
only a few percent of accelerated particles (which contains only a 
few percent of the infall energy) is injected into the intergalactic 
space. Here we must note that the escape process cannot be energy 
dependent, because then the spectrum could not match the observa- 
tions. In case of galactic cosmic rays which travel against the solar 
wind to reach the earth, strong attenuation is severely energy de- 
pendent, and so the solar wind does not provide a good example 
apparently. 

For clusters which are connected with filaments or super- 
galactic planes, the particles can escape from the clusters along 
the filaments/sheets since the field lines are mostly likely parallel 
along those structures according to the flow velocity inside the fila- 
ments/sheets (see, for example, Fig. 1 of KRJ96). In such a picture 
the low apparent efficiency is simply a measure of the difficulty 
of feeding the particles into the sheets from the surrounding of a 
cluster. If this conjecture is right, the arrival direction of UHECRs 
should show a correlation with the Supergalactic plane, and so iden- 
tifying a source cluster with a particular CR event may become less 
stringent. The particles can also leak when the accretion shock be- 
comes unstable and breaks up temporarily due to strong activities 
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Figure 5. Estimated proton flux from a cosmological ensemble of clus- 
ter accretion shocks. The top panel shows the models with Jokipii 
diffusion. The solid line is for a = 4.0, e = 5x 10~ 5 , dotted line for 
a = 4.1, e = 2xlCr 4 dashed line for a = 4.2, e = 10~ 3 and dot-dashed 
line for a = 4.4, e = 10" 1 . The bottom panel shows the results for Bohm 
diffusion. The solid line is for a = 4.0, e = 10~ 4 , and dotted line for 
a = 4.3, e = 10~ 3 . In all models we set f B = 1, m = 0, Cl b = 0.06 
and h = 0.75. The data points represent the all-experiment data set col- 
lected by T. Stanev (squares, see Rachen, Stanev & Biermann 1993), and 
the light component inferred from the Fly's Eye composition measurements 
(triangles, ibid). 



Figure 6. Estimated proton flux from a cosmological ensemble of cluster 
accretion shocks. For all models here Jokipii diffusion is assumed and e = 
5xl0~ 5 , a = 4.0, Q h = 0.06 and h = 0.75. The top panel shows the 
models with fg = 1, but with different values of m = —1,0, +1 (dotted, 
solid, and dashed lines, respectively). The bottom panel shows the models 
with m = 0, but with different values of fs = 0.1, 0.5, 1.0 (dashed, 
dotted, and solid lines respectively). The data points are the same as in Fig. 
5. 
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of AGNs inside the cluster. This could be the only possible escape 
route for the particles accelerated by isolated clusters which are not 
part of superclusters. The issue is closely related with the topol- 
ogy and the turbulent nature of the magnetic fields in the sheets, 
filaments and around clusters and will be discussed in upcoming 
papers. 

3.2 The magnetic fields outside clusters 

For our model to work, the magnetic field just outside clusters of 
galaxies must be of order 0.1 microgauss, but cannot be less, since 
then the Larmor radius condition would fail, i.e. the particles could 
not be contained. Does this imply that we require a primordial mag- 
netic field? 

We note that in the evolution of clusters, repeated formation of 
giant radio galaxies is likely, and they can break out and pollute the 
environment over very large distances with magnetic fields. There- 
fore, with such a model tracing the large scale magnetic fields to 
very large radio galaxies, one would not need a primordial mag- 
netic field. 

Of course, a primordial magnetic field would have to be 
weaker than the observational constraints, given by the cosmologi- 
cal limits on the Faraday rotation of distant radio sources; a primor- 
dial magnetic field is likely nearly homogeneous inside the bub- 
ble structure of the galaxy distribution, and so an upper limit of 
~ 200 picogauss would obtain, with a substantial strengthening to- 
wards the spatial environment near to clusters. If we use an inverse 
Parker wind model (as an analogy to the expansion of the gas from 
a radio galaxy) to do such an extrapolation, i.e. scale the magnetic 
field strength inversely with the distance to a cluster, then such an 
upper limit is just about compatible with the magnetic fields re- 
quired by our model. 

In conclusion, we cannot decide between the two possibilities, 
a primordial magnetic field, and one that derives from radio galax- 
ies, or some other sources such as normal galaxies, even on the 
large scale. 

3.3 Cosmic rays inside clusters 

The accretion shock produces cosmic rays scaling with the overall 
energetics of the cluster, but with only some fraction of the shock 
kinetic energy; thus the energy density of cosmic rays inside clus- 
ters is expected to be less than 10% of the thermal energy content in 
our model of the intracluster medium. On the other hand, the cos- 
mic rays injected from radio galaxies may produce magnetic field 
strengths and energy densities of cosmic rays close to equipartition: 
EnBlin et al. (1996) proposed that radio galaxies can push the cos- 
mic ray content in clusters to the stability limit such as is believed 
to happen in the interstellar medium due to supernova explosions 
(see, e.g. Parker 1969 for a review). As demonstrated by EnBlin et 
al. (1996) these two alternatives may be decidable in the near future 
with gamma-ray observations due to the pion decay production in 
pp-collisions. 

3.4 Arrival directions of UHECRs 

There is a study of the Haverah Park events by Stanev et al. (1995), 
which shows that the arrival directions of the UHECRs are not uni- 
form, but seem to show positive correlations with the supergalac- 
tic plane; a related study of the AGASA events by Hayashida et 
al. (1996) supports a connection between the supergalactic plane 



and the arrival directions of at least a significant fraction of UHE 
CRs. This supports models for the extragalactic origin of the high- 
est energy cosmic rays, especially the ones in which the sources are 
likely to be associated with the supergalactic plane. They include 
radio galaxies and cluster accretion shocks. We have looked at the 
distribution of local clusters (z < 0.3) and generated a map of the 
cluster accretion shocks weighted by the particle flux above a few 
times 10 EeV (Kang, Rachen & Biermann 1996). The map shows 
the similar degree of correlations with the supergalactic plane as the 
arrival directions of the UHECRs found in Stanev et al. (1995) and 
Hayashida et al. (1996). This is consistent with the fact that rich 
clusters beyond the local supergalactic plane (d > 30ft _1 Mpc) 
such as the Coma (8.3keV), Perseus (6.3keV), 3C129 (6.2keV), 
A3571 (7.6keV), and Centaurus (3.9 keV) clusters are at low su- 
pergalactic latitudes. These rich clusters are dominant contributors 
of UHE CRs according to our model. The Virgo cluster is the most 
prominent cluster in the Local Supercluster, but its temperature is 
only only 2.4 keV. Thus its contribution through the accretion shock 
is important only for energies below GZK cutoff; of course, it may 
have an additional contribution from its radio galaxies. 



3.5 The most energetic events 

Concerning the events above 100 EeV, the present model can give 
a straightforward explanation for the enhancement of highest en- 
ergy events detected by the Haverah Park experiment at high north 
galactic latitudes (see Stanev et al. 1995, and references therein), 
since the Virgo and Coma clusters are in that direction. On the 
other hand, the applicability of the hot spot acceleration model to 
the center-brightened radio galaxies M87 harbored by the Virgo 
cluster is not yet clear. Whether the 320 EeV Fly's Eye event may 
be explained in the context of cluster accretion shock accelera- 
tion, remains to be investigated. The AWM7 cluster (T = 4.0 keV, 
z = 0.0176) might be the closest rich cluster to the general direc- 
tion of this event. Alternatively, there is a candidate radio galaxy, 
3C134; its redshift is unknown due to obscuration, but its distance 
can be estimated from radio size-luminosity relations to be in the 
range 30— 300 Mpc, thus it may well be the closest FR-II radio 
galaxy (Rachen 1995). This estimate has been recently confirmed 
by an infrared detection of the central galaxy (Hartmann 1996). 
All in all, since the maximum energy of acceleration may vary 
strongly between individual objects in both models, one may have 
to consider both possibilities in the explanation of the highest en- 
ergy events. 



4 CONCLUSIONS AND OUTLOOK 

The studies on the expected CR spectrum and distribution from 
galaxy cluster accretion shocks seem to give strong evidences for a 
significant UHECR contribution, which could reach up to the high- 
est observed energies. Some aspects of the model need further de- 
tailed considerations, however, since it has to rely on several some- 
what speculative assumptions. A first question could be if one can 
prove observationally the existence of large-scale accretion shocks 
around clusters. Although the shocks and the hot component of 
IGM can be clearly seen to exist in most of cosmological hydro- 
dynamic simulations, independent of the details of cosmological 
models; it is observationally challenging to detect such shocks be- 
cause the hot gas does not radiate or absorb photons at a detectable 
level. But the hot gas of T = 10 5 -10 6 K heated by the large-scale 
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accretion shocks seems to be a major component of the intergalac- 
tic medium in the simulated universes (Nath & Biermann 1993; Os- 
triker & Cen 1996). One can also ask if the model would work for 
lower density universes such as open CDM models and flat CDM 
models with cosmological constant, because it is expected that the 
accretion flows are weaker in such cosmologies. 

Another crucial assumption for our model is the magnetic field 
in equipartition with the thermal energy in the vicinity of accretion 
shocks. The general topology and turbulent nature of the large-scale 
magnetic field are also closely related with our particle diffusion 
model and the escape process. We note that the Jokipii expres- 
sion for the possible range of the diffusion coefficient in an oblique 
shock geometry has been derived going to the absolute limit of what 
is geometrically possible; it may be an extreme overstatement of the 
limit relevant in nature. Therefore, we emphasize that we use this 
as an assumption going into the modelling. On the other hand, con- 
sidering the loss time scale and comparing the two limits, Bohm 
and Jokipii (Fig. 2), the maximum energy possible is only moder- 
ately changed, even when we are 1/3 of the way between Jokipii 
and Bohm limits. Therefore, we conclude that the speculative na- 
ture of the Jokipii limit has relatively little influence on our essential 
result, as long as nature allows us to go two thirds of the way to- 
wards Jokipii's limit. This demonstrates that the theoretical studies 
on how a magnetic field is generated and amplified during the struc- 
ture and galaxy formation (e.g. Kulsrud et al. 1996; Ryu, Kang & 
Biermann 1996; Biermann 1996), as well as the observational stud- 
ies on how the real magnetic field in the IGM is distributed are very 
important to further development of our model. 
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